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We present a thorough theoretical study of ordering phenomena in nitride ternary alloys Ga^Ini-^N, 
£SJ ALrlni-ajN, and ALjGai-^N. Using the Monte Carlo approach and energetics based on the Keating 

y—i model we analyze the influence of various factors on ordering in bulk crystals and epitaxial layers. We 

characterize the degree of both short-range order (SRO) and long ranger order (LRO) for different 

compositions, temperatures and for substrates associated with different epitaxial strain. For the 
^ description of the SRO, the Warren-Cowley parameters related to the first four coordination shells 

are used. The LRO is detected by means of the introduced sim-LRO parameter, based on the Bragg- 
Pj^ Williams approach. The description of the observed long-range ordering patterns and conditions for 

their occurrence follows. 
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I. INTRODUCTION 

Ordering in semiconducting nitride alloys is recently 
an intensively investigated issue. Both experimental as 
well as theoretical activities deal with mixtures of A1N, 
GaN, and InN, which are important from the applica- 
tion perspectives such as optoelectronics, high-power and 
high frequency electronics, and sensors. There are three 
essential issues that trigger such high research interest. 

First, there is much experimental controversy over 
which ordering patterns can be obtained in nitrides and 
under what conditions do they occur. Actually, var- 
ious modes of ordering are reported to be observed 
in the nitri de s amples. These include clus tering and 
precipitation LliSJ compositional modulation J^S or uniform 
random alloy possibly with some degree of short-range 
ordering!^! However, many of the aspects here still 
remain unclear. For example, recently Galtrey and 
coworkers^ on the basis of three dimensional atomic 
probe (3DAP) measurements, concluded that there is 
no evidence of clustering in Ga^In^^N except for the 
natural spatial fluctuations of composition. Based on a 
previous report by SmeetonP they suggested that ob- 
served "clusters" could be artifacts related to the radia- 
tive damage of the sample during high electron resolution 
transmission microscopy (HRTEM) observation. How- 
ever, this argumentation met w ith serious critique from 
TEM experimental groups) 9 ^ 10 ! pointing out that the ra- 
diative damage impact can be highly minimized through- 
out the experiment and that during such careful measure- 
ments clustering of In atoms still remains visible. In the 
presence of such ambiguities, modeling could provide a 
valuable insight into the nature of ordering occurring in 
nitrides. 

The second reason is related to the fact that ordering 
influences the electronic structure and the optical proper- 



ties of alloys. These properties, including band gap, car- 
rier localization, mobility, etc., are in turn crucial from 
the applications viewpoint, which include laser diodes 
(LDs), light-emitting diodes (LEDs), high-electron mo- 
bility transistors (HEMTs), sensors, etc. It turns out 
that even local short-range ordering significantly influ- 
ences optical properties of nitrides as well as other semi- 
conducting systems, see e.g. the work of Bellaiche et 
aZfEQ Of course, global long-range ordering has also a 
pronounced impact on optical properties in nitride al- 
loys, as it was demonstrated, e.g., in Ref. [HJ Recently, 
also Gorczyca and coworkers published a series of papers 
comparing electronic structure of small supercells calcu- 
lated assuming either clustered or uniform distribution 
of cationsPSffi^ Their studies show significant differences 
for these extreme distribution cases, for materials includ- 
ing ternary nitride alloys Ga^Ini^N, Al^Ini^^N, and 
AljJni-^N as well as selected compositions of quaternar- 
ies AL J Ga !; In 1 _ :I ._j / N. Yet, another important example of 
ordering significance for electronic properties is related 
to the luminescence intensity of In containing nitride 
samples. There are indications that excitonic recombi- 
nation occurring at In microclusters cou ld significantly 
contribute to the luminescence signalP^E^ Therefore, the 
presence or absence of In clusters should have consider- 
able impact on efficiency of light-emitting devices. 

The third issue has methodological origin. The type 
and degree of homogeneity in the considered alloy influ- 
ences the range of methods that can be applied to mod- 
eling its properties. Many theoretical approaches assume 
the perfectly random uncorrelated alloy, characterized by 
both short- and long-range order parameters equal to 
zero. Examples of such methods include virtual crystal 
approximation (VC A) , coherent potential approximation 
(CPA), special quasirandom structures approach (SQS), 
etc. The question about the presence of correlations and 
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ordering within nitride alloys is, therefore, also very vital 
in terms of accuracy and applicability of various mod- 
eling schemes. Moreover, the detailed knowledge about 
the ordering patterns could provide more realistic atom- 
istic configuration for input to semiempirical electronic 
structure computation methods such as tight-binding or 
empirical pseudopotential schemes. Thus, the overall re- 
liability of nitride alloys modeling can be seriously im- 
proved, once more detailed knowledge about ordering is 
available. 

Keeping in mind the above reasons, we try to shed 
some light on the ordering phenomena occurring in ni- 
trides. We carry out the detailed studies of ordering 
in different nitride ternary alloys including Ga^Ini-^N, 
AlaJni-^N, and A^Gai^N. We verify how various fac- 
tors could influence the ordering in these systems. The 
crucial variables include composition, temperature, and 
epitaxial strain. Our simulations are performed under 
assumption of the thermodynamics equilibrium. There- 
fore, we do not directly address the effects specific to the 
growth method. The growth of nitride layers is often per- 
formed using methods operating in non-equilibrium con- 
ditions such as molecular beam epitaxy (MBE). More- 
over, usually this epitaxial growth takes place on the 
surface. Therefore, the features related to surface or- 
dering and details of growth process could in principle 
remain "quenched" in the sample, forming a metastable 
state. Phenomena of this kind are not directly included 
in our modeling. The main computational tool of our 
investigation is the static Monte Carlo method. This 
is not the only option. There are also kinetic studies 
available in the literature for nitrides.^ The energet- 
ics in our calculations is computed on the basis of the 
Keating valence force field model and, for the sake of 
simplicity, we focus on the zinc-blende structures. The 
studies are carried out in the lattice-coherent thermo- 
dynamics regime, i.e., we assume that the alloy con- 
stituents form common lattice and can not decay into 
separate components, each with its own independent lat- 
tice parameter. This assumption was often not included 
and there were numerous studies examining the behav- 
ior in the so called lattice-incoherent case. These ap- 
proaches rely mainly on the free-energy difference be- 
tween alloy and binary constituents relaxed to their own 
lattice constants. Examples of this type of cal culation for 
ternarjffSHSU anc j f or quaternary nitrides^l22l are present 

in the literature. These studies provided qualitatively 
similar phase diagrams with unimodal curves separating 
the uniform region in the composition-temperature pa- 
rameter space. Moreover, they predicted very high criti- 
cal temperatures Tm, above which the alloy components 
are fully miscible. For the most analyzed case of the 
ternary Ga x liii_ x N, typical reported values of Tm are in 
the range of T M w 1500-2000 K (see the work of Liu 
and Zunger 3 -^ for a thorough review of different model 
findings) . These high Tm values yield very low estimates 
for the maximum In composition, which can be incorpo- 
rated in Ga x Ini_2;N in typical growth temperatures, be- 



fore phase separation occurs. On the other hand, there 
are experimental findings showing that it is possible to 
grow Ga x lni_ x N samples containing up to 20-30% of In 
without triggering significant clustering, which is in dis- 
agreement with aforementioned predictions of the lattice- 
incoherent case. Therefore, recently also a different ap- 
proach was introduced EDM] where the lattice coherence 
is assumed. It turned out that in this case Tm decreases 
to values below typical growth temperatures. Instead of 
the phase separation, the random alloy phase with some 
degree of short-range order is predicted. Therefore, in 
our study, we focus on lattice-coherent case, which seems 
to correspond better to the alloys obtained during epi- 
taxial growth. Moreover, we also analyze the influence 
of the stress related to different epitaxial substrates em- 
ployed, as this can be an impor tant factor suppressing or 
triggering rdering P 2 | 23 | 28 | 31 | 3 S 

In order to characterize the degree of ordering we ex- 
amine equilibrium structures obtained from our Monte 
Carlo simulations. To quantify the degree of short-range 
order in the generated samples we use the Warren-Cowley 
family of parameters. For description and detection of 
long-range ordering we employ the approach based on the 
Bragg- Williams long range order measure. We present 
a complete and systematic study of these parameters in 
the whole composition range for all ternary combinations 
Ga^Ini-zN, ALTni-^N, and Al^Gai-^N. We also exam- 
ine different temperatures, beginning from growth tem- 
perature range to higher values which could correspond 
to the sample annealing conditions. Moreover, we exam- 
ine how the above ordering metrics are influenced by bi- 
axial strain related to application of a different substrate 
for growth process. 

The paper is organized as follows. We start with a 
short overview of the short- and long-range ordering con- 
cepts in Sec.[n] The commonly employed Warren-Cowley 
and Bragg- Williams parameters are also briefly discussed 
there. Next, in Sec. |III[ the details of the employed 
computational approach are explained. Section |IV| de- 
scribes how the order parameters are computed from the 
Monte Carlo data. In particular, it introduces the sim- 
LRO parameter, constructed on the basis of the Bragg- 
Williams characteristic. The sim-LRO metric proposed 
in this work is a handy indicator of the long-range order- 
ing emergence during Monte Carlo simulations. The re- 
sults for the ternary bulk alloys Ga^Ini-^N, ALTni-^N, 
and ALjGai-zN are described in Sec.|Vj The ordering in 
biaxially strained epitaxial Ga^Ini^^N layers is studied 
in Sec. |VI| Finally the paper is concluded in Sec. |VII| 



II. QUANTITATIVE DESCRIPTION OF 
ORDERING 

The concepts of order and disorder are intuitively eas- 
ily understood. However, there exist many ways of con- 
verting this intuition to numbers describing the degree 
of ordering in crystals. Therefore, to provide a neces- 
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sary background and fix the notation, a brief description 
of the concepts used throughout the forthcoming anal- 
ysis of ordering in nitrides is given in this section. For 
more thorough insight into the ordering phenomena, see 
classical monographs of ZimanpS! DucastelleJ^ or a very 
informative introductory paper by Klein.^1 



A. Long-range order vs short-range order 

When analyzing the ordering in alloys, it is important 
to distinguish between two different aspects of this phe- 
nomena, namely, the short-range order (SRO) and the 
long-range order (LRO). The term short-range order is 
used to describe the preference of certain types of atoms 
to reside near each other. This effect manifests itself in 
the form of statistical correlations between occurrences 
of atomic types on their respective coordination shells. 
Usually, alloys exhibit non-zero degree of short-range or- 
der, as a result of energetic preference towards particular 
atomic arrangements. The long-range order is related to 
the development of a global pattern spread throughout 
the whole crystal. It is worth underlining that in the tem- 
peratures above K the observed pattern is never perfect 
and contains some misplaced atoms. Both types of or- 
dering can be described in a quantitative manner. In the 
following we recall the definitions of the commonly used 
Warren-Cowley SRO parameter and the Bragg- Williams 
LRO characteristic. 



B. Warren-Cowley short-range order parameter 

One of the most commonly used measures of SRO is 
the Warren-Cowley family of parameters rPSl In multi- 
component alloys, these SRO parameters between com- 
ponents A and B are defined on the basis of conditional 
probability of finding an A-type atom in the specified 
coordination shell of a B-type atom: 
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P (A atom in shell i of site j \ B atom on site j) 
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(1) 



where ca denotes the concentration of A-type atoms in 
the system. It is straightforward to note that the relation 
Tab < indicates the preference to AB neighborhood 
on the ith coordination shell, whereas T A g > indicates 
anti-preference. For the ideal uncorrelated random alloy, 
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and that for each of the components P the following sum 
rule holds: 



B 



cbT® b = 0. 



(3) 



Therefore, for the alloy Aa;Bi_ x , we can in principle 
define four different SRO measures Taa, Tab, Tea, 
and Tbb- However, because of one symmetry relation 
[Eq. ([2|] and two sum rules [Eq. ^3|], actually only one 
SRO parameter is really independent. In this paper, we 
use standard convention and focus on Tab- It is worth 
underlining that, even though technically we deal with 
ternary alloys A^Bi^N, they can be described identi- 
cally as in the binary case. This is because nitrogen 
lattice remains mostly unaffected, as the probability of 
substitution in the nitrogen site is much lower than ex- 
changes in cationic lattice because of the energetic rea- 
sons. Therefore, we focus only on the cationic fee sub- 
lattice and investigate in a great detail order parameters 
related to the first four coordination shells within this 
lattice. 



C. Bragg- Williams long-range order parameter 

The long-range order must be specified with respect to 
a certain pattern. To calculate it, e.g., for A 2 .B 1 _ 2 . binary 
alloy we have to know which sites should be occupied by 
atoms of type A, type B, etc. Once we know the spatial 
pattern (PAT), we can define the Bragg- Williams long- 
range order parameter as 



,(PAT) 
j(PAT) _ Ja 



CA 



1 - C A 



(4) 



For the SRO parameter, a key role was played by the 
conditional probability, whereas for the LRO measure 

the most important variable is / A PAT ' ) . It denotes the 
fraction of A sites from pattern PAT, which are actually 
occupied by A type atoms in the structure under consid- 
eration. Symbol ca stands here for concentration of A 

{ PAT) 

atoms in the sample. If 5 A ~ 1, it indicates that the 
location of A atoms in the structure is similar to pattern 

(PAT) 

PAT. If S A = 0> ^ 

means that no significant similar- 
ity to PAT in terms of A-atom locations was detected. 
From sf AT) < 0, it can be deduced that A atoms in 
the structure avoid A sites from PAT. For the binary al- 
loys A x Bi_ x , the 5 A PAT ' 1 = S'b PAT ' 1 , so we have only one 
independent Bragg- Williams LRO parameter. For the 
technical reasons, the direct use of the Bragg- Williams 
parameter for the quantification of the LRO in Monte 
Carlo simulation data is difficult. Therefore, to detect 
the presence of the LRO in our results, we develop a de- 
rived quantity, the sim-LRO parameter. It is described 



in greater detail in Sec. IV 



III. COMPUTATIONAL DETAILS 

All calculations, were performed for the zinc blende 
6x6x6 supercell containing 1728 atoms. As resulted 
from our tests, this size of the supercell ensured good 
convergence of the examined order parameters. See 
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Convergence of r' 1 ' for Ga 25 ln 75 N in MxMxM supercell 
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FIG. 1. (Color online) Sample convergence test — the de- 
pendence of the Warren-Cowley short-range order parameter 
^gLiu on tne M x M x M supercell size in bulk Gao.25Ino.75N. 
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FIG. 2. (Color online) Differences between short-range order 
parameters estimation from simulations powered by two dif- 
ferent random number generator — RanLux from the GSL 
library and Wichmann-Hill from the ACML package. Simu- 
lations were carried out for the Ga^Ini-ajN on InN substrate 
in T = 873 K. 



Fig. [T] for a sample test, showing the dependence of 

(i) 

-[gain in bulk Gao.25Ino.75N on the simulation supercell 
size. When it comes to the applied simulation tech- 
niques, the standard static Metropolis algorithm within 
the NVT ensemble was applied. The concentration of 
species within alloy was held constant during the sim- 
ulation, atomic shifts and cationic exchanges were em- 
ployed as trial moves. Because of the employment of 
atomic shifts the lattice vibrational effects are directly 
included in the performed calculations. The magnitude 
of atomic shifts was adjusted so that exchange proba- 
bility was equal to 0.5, as reco mmen ded for the most 
efficient phase space exploration! 39 * 40 ! The exchange tri- 
als were performed five times more often than atomic 
shifts, since in the majority of the cases their probabil- 
ity was much lower than for shifts. The total length of 
simulation was 8.3 million Monte Carlo sweeps per con- 
centration. For each run, we allowed 0.4 million Monte 
Carlo sweeps for equilibration, before gathering the sim- 
ulation statistics. The energy calculations were carried 
out using Keating modeP^ and its parametrization for 
nitrides described in the earlier work of the authors.^ 
Periodic boundary conditions were used throughout all 
presented simulations. For every type of ternary alloy, 
A^Bi-rcN a series of simulations with 17 different con- 
centrations x — 0.056, 0.111, 0.167, . . . , 0.944 was per- 
formed. For the initial conditions, random and uncor- 
rected structures were generated with lattice constant 
optimized to their energy minimum. For the pseudoran- 
dom number generator, the RanLux method was used as 
implemented in the GSL library.^! Also tests with dif- 
ferent algorithms were performed, showing that the sim- 
ulation results are insensitive to pseudorandom number 
generation technique employed. Figure [2] displays the 

difference between the parameters computed for 

Ga^Ini-^N on InN using the aforementioned RanLux al- 
gorithm and Wichmann-Hill method implemented in the 
ACML librarypS This difference was calculated according 



to the formula 

DT<£i In = T« In (RanLux) - rg^ In (Wichmann-Hill) . (5) 

The observed values of DT^^ are small. As it can be 
seen from Fig. [2] they do not exceed a few percent. 



IV. DETERMINATIONS OF ORDER 
PARAMETERS FROM SIMULATION DATA 

(i) 

The computation of Warren-Cowley parameters T^g 
from structures obtained in Monte Carlo simulations is 
relatively easy. It is enough to estimate the conditional 
probability from formula ([I]) by average percentage of 
type A atoms present in the ith shell of the B atoms in 
the set of structures generated during each Monte Carlo 
run. However, as opposed to the short-range order, the 
situation for the long-range ordering is not so straight- 
forward. To be able to estimate the Bragg- Williams 
measure «S'^ PAT ' ) , the pattern PAT has to be known in 
advance. Its correct selection is crucial for the detec- 
tion of the long-range order. When calculating S A 
for a highly ordered structure, however, very different 
from assumed PAT, one may easily obtain values close 
to zero, indicating the absence of LRO. Moreover, even 
if the exact pattern is known, one should check for all 
its symmetry-equivalent variations, since the value of 
S^rat) clearly depends on assumed orientation of PAT. 
Therefore, in principle, all variants of PAT produced by 
symmetry operations should be checked, as the variant 
developed during simulation is usually a result of random 
fluctuations and can not be predicted a priori. However, 
it would be very useful to have a quick way of checking 
whether long-range order of some type develops during 
Monte Carlo simulation or not, preferably without the a 
priori knowledge about ordering pattern. 
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It turns out that, in order to detect the presence of 



(PAT) 
A 



long-range ordering, quantities derived from 5 
more useful. Instead of a priori selecting ordered config- 
uration PAT, it is helpful to calculate the Bragg- Williams 
parameter with one of the structures obtained in the de- 
veloped stage of simulation used as pattern PAT. It is 
important that this reference structure PAT is selected 
far after the thermalization phase. The quantity calcu- 
lated with this method we will call the simulation LRO 
parameter (sim-LRO) and denote it as 5a- The typical 
behavior patterns of 5a during Monte Carlo simulation 
run are displayed in Fig. [3] They are taken from calcu- 
lations for Ga^Ini-^N on InN substrate, which are dis- 
cussed in greater detail in Sec. |VI| It turns out that three 
types of behavior for 5a can be distinguished. First, if 
no LRO is present then Sa simply fluctuates around 
value, corresponding to unordered alloy. This is il- 
lustrated in Fig. [3] (a). Second, if certain LRO pattern 
was developed during simulation, the successive struc- 
tures fluctuate around it generating values of 5a that 
are considerably different from 0. This is depicted in 
Fig. [3] (b). Third, sometimes different arrangements of 
atoms could repeat during simulations, e.g., related to 
different orientation or shift of LRO pattern. It will cor- 
respond to series of peaks visible in Sa values during 
simulation time, see Fig. [3] (c). These peaks correspond 
to moments when the structure is similar to the one se- 
lected for comparison. Even though, the average values 
of 5a throughout the whole simulation could be in this 
case close to 0, this behavior also indicates the develop- 
ment of LRO or precipitation. Therefore, to detect this 
type of ordering, except the average value of 5a we also 
analyze the spread of observed 5a values. Convenient 
measure here is the difference of percentiles q related to 5 
and 95 percent, estimated on the basis of 5a population 
generated during Monte Carlo (MC) sampling 



(a) 



S Ga forGa 017 ln 83 N on InN 



A5a = go.95 (5a) — 90.05 (5a)- 



(6) 



This quantity provides the interval width which includes 
90% of the observed data. The occurrence of high values 
for either average sim-LRO parameter 5a or its spread 
A 5a, enables for convenient detection of long range or- 
dering in MC simulation results. The average value of 5a 
together with the percentiles forming A5a are marked 
with dashed lines in Figs. [3] (a)-(c). 



V. ORDERING IN THE TERNARY BULK 
ALLOYS 

In this section, the ordering in the bulk nitride ternary 
materials Ga x Ini_ x N, Al^Ini-^N, and Al^Gai-^N is ex- 
amined in a great detail. In the analyzed cases, no long- 
range ordering was observed. The sim-LRO parameter 
for all presented cases fluctuated around zero, similarly to 
the situation displayed in Fig. [3] (a). Therefore, we focus 
here on the short-range ordering. The Warren-Cowley 
SRO parameters for the first four coordination shells in 
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FIG. 3. (Color online) Different types of behavior of Sea 
during simulations for for Ga^Ini-^N alloy on InN substrate. 
Three gallium concentrations x are presented: x = 0.17 (a), 
x — 0.50 (b), x — 0.72 (c). Dashed lines correspond to the 
average value of Sgs, and its spread A5G a . 



the cationic lattice are examined in the whole concentra- 
tion range. Moreover, this study covers a few different 
temperatures, beginning roughly from the lower end of 
growth temperatures range, i.e., 873 K, and including 
higher values up to 1673 K. This upper limit could corre- 
spond to annealing during the post-processing phase. In 
the following the results for each of the ternary combina- 
tions, namely, Ga^Ini^N, Al x Ini_ x N, and Al^Gai-^N, 
are described. 



A. Ordering in Ga^Ini-^N 

Ga^Ini-^N is up to now the most investigated of all 
ternary combinations of nitrides. Currently, it is the 
main ingredient of the active region in blue-green op- 
toelectronic devices. The dependencies of r^ In , r^ In , 

^Gain' an d TQ^in on composition for different tempera- 
tures are presented in Fig. [4] The most interesting com- 
position range, corresponding to violet-blue-green wave- 
length 400-570 nm is gray shaded on all graphs. It is 
easy to observe that r^ In and T^l ln are negative in 
the whole concentration range and for all examined tem- 
peratures. This indicates the preference toward Ga-In 
neighboring on the first and fourth coordination shells. 
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(2) (3) 

Conversely, and remain positive, which in- 

dicates that Ga atom is more likely to have another Ga 
atom on its second and third coordination shell, than it 
would result from concentration x. This sign pattern is 
in qualitative agreement with recent results of Chan and 
coworkers^ obtained on the basis of the cluster expan- 
sion model. This —/ + / + /— sequence observed for the 
rW parameters is usually interpreted as manifestation 
of stability for the chalcopyr ite st ructure, which exhibits 
the same SRO signs pattern 

For each dependence of the short-range order param- 
eter on gallium concentration x, the third order polyno- 
mial was fitted. This type of curve is selected as the 
simplest model capable of describing observed shapes 
with asymmetric maximum/minimum. These fits are 
presented as continuous lines in Fig. [4] It turns out that 
all short-range order parameters dependencies have sin- 
gle extremum, corresponding to maximum absolute value 
of r parameter. On the basis of obtained polynomial fits, 
concentrations x^ t corresponding to this maximum or- 
dering are determined. These extrema are marked with 
open symbols in Fig. [4] The numerical values of ob- 
tained Xcxt are presented in Table [i] For the nearest 
neighbors and next nearest neighbors, the cc^t is located 
at lower gallium concentration roughly around 35%, or in 
other words at high indium content 65%. These composi- 
tions are outside of the most interesting violet-blue-green 
range. For the third and fourth coordination shells, these 
extrema move toward approximately 50% compositions. 
Also the overall curve shape gets more symmetric around 
extremum for the case of r^L^ and r^ In than it is 
in the case of the first and second shell. Nevertheless, 
these higher-order parameters should have weaker influ- 
ence, e.g., on electronic properties of Ga x Ini_a;N than 



.(i) 



(2) 

GaIn and rQ^ In . Therefore, we conclude that the re- 
gion, where effects related to ordering are the most pro- 
nounced is located around x ~ 0.35. It is also worth 
noting that the position of x^ t is virtually independent 
on temperature, and the difference 
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(T = 1673K) 



(T = 873K) (7) 



for all coordination shells i = 1, ... 4 does not exceed 0.04. 
We also calculate the difference of extremum values for 
each r« 



AT. 



r« (T = 1673X) - TZ(T = 873K). (8) 



(i) 



This quantity indicates that the ordering decreases 
roughly by half when the temperature changes from 
T = 873 to T = 1673 K. All the above findings are sum- 
marized in Table |U 

As far as modeling of electronic structure for 
Ga^Ini-^N is concerned, the above results indicate that 
it would be the most interesting to focus on structures 
with having successive signs —/ + / + /— as an 

input to density functional theory, tight-binding or em- 
pirical pseudopotential schemes. Methods which assume 
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FIG. 4. (Color online) The short-range order parameters in 
Ga^Ini-ajN as a function of composition for different tem- 
peratures. Note that the scale on all graphs is the same to 
allow for direct comparison of ordering magnitude in differ- 
ent coordination shells. Continuous lines represent third order 
polynomial fits to simulated data. The extrema for presented 
fits are marked with open symbols. The composition range 
approximately corresponding to violet-blue-green wavelength 
of emitted light (400-570 nm) is gray shaded. 
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^Gain = such as, e.g., special quasirandom structures 
(SQS) or virtual crystal approximation (VCA) could lead 
to systematic inaccuracy, caused by neglecting of the 
SRO. The range, where these inaccuracies should be the 
most pronounced, is determined by extrema for the low- 
est coordination shells SRO parameter, namely around 
indium concentration of 65% or so. On the other hand, 
when considering alloys corresponding to the violet end 
of spectra, say around 10% of indium, the influence of 
ordering effects should be much less pronounced. The 
comparison of the above results with experiment is not 
straightforward. There is a heated debate about the 
range of concentrations and, more generally, conditions 
that lead to In clustering in Ga^In^^N samples. This 
issue so far is by no means clarified, even on the phe- 
nomenological level. Our computations, to a certain ex- 
tent, support these findings that report no clustering 
or long-range ordering in the samples even with high 
In concentrationPtH I n our simulations we show that 
under the lattice-coherent assumption, thermodynamics 
does not prohibit obtaining mixtures without precipi- 
tates or LRO in very broad indium concentration range 
and in temperatures higher or in the region of the typ- 
ical growth conditions. However, one has to stress the 
fact, that this state could be difficult to obtain exper- 
imentally, mostly due to artifacts related to epitaxial 
growth methods. These methods, such as MBE or metal- 
organic vapor phase epitaxy (MOVPE), take place on 
the surface, in conditions at least partially correspond- 
ing to non-equilibrium. Therefore, some features specific 
to growth process could form metastable state and re- 
main "quenched" in the samples. To fully describe these 
phenomena, the detailed growth models should be devel- 
oped. This is, however, far more difficult than already 
not-easy computation of equilibrium properties. Yet an- 
other issue not addressed on this stage of modeling is 
the presence of strain, both local and global. By local 
strain we mean deformations related to defects, in partic- 
ular dislocations, whereas global strain occurs in thin lay- 



TABLE I. The summary of extrema for T(x) dependencies in 
GaaJni-zN presented in Fig. [4] Concentrations of maximum 
ordering x^ t , together with corresponding extremum value 
of at the lowest examined temperature T = 873 K are 

provided. The influence of temperature is illustrated using 
the difference Arr^t between extremum concentrations at T = 



1673 K and T = 873 K [see Eq. 
SRO parameters Ar^t [see Eq. 



7 1] and extremum values of 
81]. 



,(0 



r«(T=873 K) Aa£&(T=873 K) Ar 



-i(i) 

■ ext 



r&L 0.35 



r (2) 
1 Gain 

p(3) 
L Gain 



0.36 
0.47 



r&„ 0.52 



-0.088 
0.053 
0.033 

-0.050 



-0.006 
0.023 
0.040 

-0.015 



0.040 
-0.025 
-0.019 

0.022 



ers and is related to lattice mismatch between substrate 
and alloy material. We reference it as global, since it is 
present in the whole volume of the considered layer. The 
influence of dislocations and other sources of local strain 
is not included in our modeling, since we assume perfect 
crystalline lattice. Taking into account such extended 
defects would require calculations with much larger su- 
percells with complicated geometry, prohibitive for the 
detailed composition and temperature scans provided in 
this study. Nevertheless, we can gain certain insight into 
the influence of strain on ordering by modeling global 
epitaxial deformation. It was already indicated in the 
literature that this epitaxial s train can ha ve significant 
impact on ordering phenomena! 22 * 23 * 28 * 31 * 34 ^ Therefore, its 
role is analyzed in detail in Sec. |VI| 



B. Ordering in Al^Ini-^N 

When speaking of ternary nitride alloys with applica- 
tion prospects, Ga^Ini-^N is not the only option. Re- 
cently, also ALjIni_ x N attracts considerable attention. 
It is particularly interesting, because by adjusting the In 
composition one can lattice match it to GaN substrate, 
producing in principle strain-free interfaces and hetero- 
junctions. The absence of strain has many advantages 
such as lack of electric field component related to piezo- 
electric effect, lower density of defects on the interface, 
etc. Therefore, such Al^Ini-^N layers are promising can- 
didates for applications in optoelectronic devices, e.g., for 
high quality factor microcavities or distributed Bragg re- 
flectors. Moreover, the AlInN/GaN junction could be an 
important building block for high electron mobility tran- 
sistors (HEMTs). For more detailed information about 
properties and applications of Al 2 .In 1 _ a; N, see a recent 
review paper by Butte and co-workers.^ 

In our study of ordering in the ternary ALTni^N, we 
observe qualitatively similar behavior as in the case of 
Ga^Ini-^N (compare Figs. [4] and [5|. Again, the same 
alternating pattern of signs for order parameters is ob- 
served, i.e., rjyj n and Ijyj n are negative, whereas r^j n 

(3) 

and r^ 1In are positive. Similarly to the previous sec- 
tion, the SRO parameter values and concentrations for 
extremum ordering r^ t and x^t are determined. The 
extremal values of order parameters were in the case of 
ALTii 1 _ ;r N larger than in the case of Ga^In^^N roughly 
by a factor of 50%, which correlates with larger lattice 
mismatch between A1N and InN than between the GaN 
and InN. For detailed comparison, see Table [TT] and the 
analogous data for Ga^Ini-^N gathered in Table [I] Im- 
portant difference in the case of ALTni-zN is that the 
composition dependencies of r^l, I^L, and r^j n 
exhibit features of bimodal shape, particularly in lower 
temperatures. Therefore, in these cases the third-order 
polynomial fits carried out in order to locate extrema, 
were performed in the narrowed range: x € [0.0,0.6] for 



.(2) 
Alln' 



and x £ [0.0,0.7] for r£/ In . The bimodal 



(4) 

Alln' 



TABLE II. The summary of extrema for T(x) dependencies 
in ALjIni-^N presented in Fig. [5] Concentrations of max- 
imum ordering a^ xt , together with corresponding extremum 
value of r^j In at the lowest examined temperature T = 873 K 

(2) (31 

are provided. Since in ALTni-^N, for r^ 1In and r A1In two 
extrema are present, both of these are quoted in the table. 
The influence of temperature is illustrated using the differ- 
ence Ax^t between extremum concentrations at T = 1673 K 



(a) 



, for Al y ln<_ x N 



and T = 873 K, [see Eq. 



and extremum values of SRO 



parameters Ar^ t [see Eq 



,(0 



r^ t (T=873 K) Aa; cxt (T=873 K) Ar 



(0 



(1) 

Alln 


0.31 


-0.131 


0.006 


0.055 


(2) 
Alln 


0.31 
0.70 


0.090 
0.067 


0.024 


-0.042 


(3) 
Alln 


0.41 
0.66 


0.058 
0.052 


0.033 


-0.034 


(4) 
Alln 


0.50 


-0.076 


0.016 


0.033 



character is particularly pronounced for r^l , 1^2 in 
the lowest temperature T = 873 K. In these cases, the 
second maxima are clearly visible and could be easily cal- 
culated from polynomial fits. The detailed information 
about all extrema is provided in Table [n) Apart from the 
bimodal character, other features observed for the SRO 
parameters in AL^Ini-aN are similar to Ga a Ini_ x N. The 
shape of extremum for 1^2 and r? lr is asymmet- 
ric and located in the high indium concentration range 
around 70%, a little bit higher than for the Ga^Ini-^N 
case. For the highest coordination shells, the extremum 
gets more symmetric and shifts toward 50% composition. 
Again for this alloy, the location of extremum concentra- 
tion weakly depends on temperature and changes maxi- 
mally around 3%, as presented in Table [5) 

Since the properties of Al x Ini_ x N lattice matched to 
GaN are of high practical importance, we gather the 
information about the short-range order parameters of 
this material in Table IIII1 In the context of wurtzite 
ALTni-^N, usually the aluminum concentration x ~ 0.82 
is quoted as corresponding to lattice matching.^ In the 
case of cubic materials and lattice parameters employed 
in this study, the Vegard's law leads to the concentra- 
tion x s» 0.79. For simplicity, we neglect here the small 
unphysical lattice bowing effect which is a feature of the 
Keating model.lSJThe composition of the Al^Ini-^N lat- 
tice matching GaN, is also marked for convenience in 
Fig. [5] It turns out that the magnitude of SRO parame- 
ters in this case does not exceed 0.05. Interestingly, it is 
the lowest for the first coordination shell, which should 
have the greatest impact on the alloy properties. When 
it comes to modeling implications, the situation here is 
analogous to Ga 2 .In 1 _ :! .N. Our calculations indicate that 
for electronic structure modeling the configuration with 
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FIG. 5. (Color online) The short-range order parameters in 
ALjJni-zN as a function of composition for different tempera- 
tures. Note that the scale on all graphs is the same to allow for 
direct comparison of ordering magnitude in different coordi- 
nation shells. Continuous lines represent third order polyno- 
mial fits to simulated data. The extrema for presented fits are 
marked with open symbols. The composition of AUIni-^N 
lattice matched to GaN is indicated with gray line. 
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SRO sign pattern —/ + / + /— should be the most inter- 
esting. The deviation from random uncorrelated alloy is 
the biggest for high indium concentration around 70%, 
therefore, in this region models assuming SRO=0 should 
have the largest systematic inaccuracy. 



TABLE III. Summary of the short-range order parameters 
for Al x Ini_a;N lattice matched to GaN, i.e., for aluminum 
concentration x ~ 0.79. 



T=873 K 



T=1673 K 



r (i) 

1 Alln 

r (2) 
1 Alln 

p(3) 
1 Alln 

r (4) 
1 Alln 



-0.014 -0.012 

0.049 0.019 

0.041 0.016 

-0.038 -0.027 



C. Ordering in AlzGai-^N 

Finally, the last possible ternary combination for the 
examined family of nitrides is Al^Ga^^N. This material 
is very well suited for high electron mobility transistors, 
operating in high-power range as well as ultraviolet light 
emitters and detectors. It is also a promising building 
block for biosensors. Since the lattice mismatch and the 
differences in force field parameters are the smallest in 
this case, also the ordering effects here are the weak- 
est. For this material, we carried out the studies only for 
two temperatures 873 and 1673 K. None of the examined 

l-AiGa va l ues exceeded 0.01. Therefore, our conclusion is 
that ALjGai-^N follows quite closely the picture of of 
uncorrelated random alloy, i.e., with both SRO=0 and 
LRO=0. 



VI. ORDERING IN TERNARY ALLOYS ON 
THE SUBSTRATE 

In the previous section we have examined ordering phe- 
nomena in bulk crystals. However, for applications in 
optoelectronic devices or sensors usually thin epitaxial 
layers grown on substrate are employed. In this part, 
we describe, how the presence of substrate associated 
with biaxial strain influences the ordering phenomena. 
We focus on the technologically most important case of 
Ga :E In 1 _ ;r N and examine its behavior on a variety of sub- 
strates. However, before moving on to the discussion of 
ordering, the basic facts about thin epitaxial layers are 
summarized below. 

The main parameter for epitaxial layer grown on the 
substrate is misht strain. For a cubic material, it is de- 
fined as 



where a su bstrate and ai ayC r denote substrate and layer al- 
loy lattice constants respectively. If we consider the sim- 
plest model, the misfit strain is compensated by purely 
elastic deformation of the thin film material. The thick 
substrate is assumed to stay undeformed in this ap- 
proach. See Fig. [6] for illustration of such a situation. 
Within this model, the deformation of layer unit cell on 
the interface (in order to match the substrate lattice) 
is compensated by relaxation in the perpendicular direc- 
tion. One can approximate the relaxation strain e re i a x as- 
sociated with accommodation of misfit strain using linear 
elasticity theory 



^rcla 



9 c 

& c misfit ■ 

Cll 



(10) 



However, for large misfit strains going beyond the ap- 
plicability range of the linear elasticity (few percent de- 
formations), this expression can be a rather rough ap- 
proximation. It is important to note that in the case 
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FIG. 6. (Color online) The epitaxial layer in the presence of 
lattice mismatch between the layer material and the substrate 
material. The case of purely elastic accommodation of strain 
is presented, i.e., without misfit dislocations (lattice-coherent 
case) . 



of lattice mismatch between the substrate and epitaxial, 
material there exists a critical thickness h, beyond which 
layer starts to relax the strain in the form of misfit dis- 
locations and other defects, instead of purely elastic de- 
formation. This means that because of the high defects 
density, epitaxial layers beyond h are virtually useless 
from the point of view of light-emitting heterostructures. 
However, estimation of the critical thickness is a difficult 
problem. Various models were proposed for this purpose. 
In this paper, in order to get a rough estimate of the criti- 
cal thickness, we employ a simple approach developed by 
Matthews and Blakeslee^ and later on employed to cubic 
nitrides.- 8 It predicts that h is given by the transcenden- 



10 



tal equation 



h = 



^substrate (1-I//4) 
4V2iT {l + v) |e misfit | 



In 



V2h 

^substrate 



(ii) 



where v is the Poisson ratio of the layer material v — 
c\%l(cx\ + C12). Note that since both misfit strain and 
Poisson ratio depend on alloy composition, the result- 
ing critical thickness also depends on the proportion of 
compounds in alloy epitaxial layer. Employed values of 
lattice constants are given in Table |IV| the composition 
dependencies of the elastic parameters in Eq. (11) were 
taken from our previous worlJ^l an d are given by 

Cll (x) = 182.23 + 105.70 x - 15.44 x (1 - x), (12) 
c u (x) = 104.78 + 25.02 a;- l.OOx(l-ar). 

For certain compositions related to large misfit strains, 
the real- valued solutions of the Eq. (11) do not exist, 



which indicates that fabricating even very thin layers 
of high crystalline quality is impossible for certain al- 
loy/substrate combinations. However, one has to bear in 
mind, that accurate modeling of the critical thickness is 
complicated task in itself. The Matthews and Blakeslee 
approach is simple, but crude model. Therefore, it gives 
more qualitative than quantitative picture of the h de- 
pendence on composition and limits of the elastic strain 
accommodation regime. 

In the present section we examine the Ga^Ini-^N 
ternary alloy grown on a large variety of possible sub- 
strates, both employed experimentally as well as possibly 
suitable in the future, including 3C-SiC, zb-ZnO, CaO, 
A1N, GaN, and InN. Currently, the most promising well 
lattice matched substrates for the zinc blende nitrides 
seem to be 3C-SiC or GaN obtained either in the form 
of thick buffer layer on a different substrate or recently 
fabricated zb-GaN free-standing crystalsP^HUl However, 
with technological progress, other options can gain sig- 
nificance in the future. The information about misfits 
bounds for each of these materials is gathered in Ta- 
bic US 

To simplify the discussion, we separate our results into 
two groups of substrates related to moderate and large 
misfit strains, since the observed phenomena in both 
cases are quite different. The first group consists of 
GaN, zb-ZnO, and CaO, whereas the second comprises 
3C-SiC, A1N, and InN. All the results have been obtained 
in T = 873 K, corresponding to a typical growth temper- 
ature range. 



A. Moderate misfits regime 

For substrates related to moderate misfit strains (i.e., 
GaN, zb-ZnO, and CaO), no long-range ordering is ob- 
served in Ga^Ini^^N, similarly to the bulk case. Since 
the intensity of ordering effects in materials is quanti- 
fied by the absolute value of short-range order parameter 



TABLE IV. Misfit strain range for various substrate applica- 
ble to GaaJni-ajN alloys. The extremum values for pure GaN 
and pure InN are provided. 



Substrate 


Lattice 
constant 


Misfit strain ranj 

x= 1 
^pure LralN J 


;e for Ga^Ini-^N 
x = 
(pure InN) 


3C-SiC 


4.360 


-3.2 % 


-12.8 % 


A1N 


4.374 


-2.9 % 


-12.5 % 


GaN 


4.503 


0.0 % 


-9.9 % 


zb-ZnO 


4.580 


1.7 % 


-8.4 % 


CaO 


4.811 


6.8 % 


-3.8 % 


InN 


5.000 


11.0 % 


0.0 % 



Gain I' 



it is useful to analyze the influence of the sub- 



strate in terms of Ar^ In defined as 

Alf aIn = ^(on substrate)| - |r« In (bulk)|. (13) 

If the presence of the substrate does not change the sign 
of 

l^Gain ( w hich is the case for examined GaN, zb-ZnO, 
and CaO), the interpretation of AFQ aIn is straightfor- 
ward. When ArQ^ In > 0, the ordering effects are in- 
creased, whereas for Ar^^ < 0, the grown thin layer 
behaves closer to the uncorrelated random alloy than the 
bulk case. 

As an example of obtained results, the case of GaN 
substrate is presented in Fig. [7| The findings for CaO 
and ZnO are qualitatively similar. The graphs illus- 
trate the comparison of the short-range order parameters 

r Gai n >_ r Gain. Ig^, and r Gai n between strained and 
unstrained cases. Order parameters are accompanied by 
the corresponding ArQ aIn values. Also, the misfit strain 
Cmisfit and the relaxation strain e re i a x dependence on com- 
position are provided for completeness. Finally, the in- 
formation about the critical thickness h calculated from 



the Matthews-Blakeslee model given by Eq. (11) is in- 
cluded. This provides insight, as to range of composition 
is of practical importance. 

The e ro i ax presented in Fig. [7] (c) has been calculated 
using two methods: first, within the theory of elasticity 
(TOE) and Eq. (10), second by relaxing the unit cell to 
the shape corresponding to minimum energy employing 
the Keating valence force field (VFF). For GaN and other 
examined substrates, both methods give similar results. 
The small differences emerged for compositions related to 
larger misfit strains. This is intuitively well understood, 
since the Keating model was constructed to recover the 
predictions of the TOE for small strains.^ We also ob- 
served that for negative misfit strains (compression in the 
substrate plane) , the theory of elasticity slightly overesti- 
mates the e re i ax , whereas for positive misfits (in the case 
of zb-ZnO and CaO substrate), the e re i ax predicted by 
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r Gain ,or < 3a x ln 1 _ x N grown on GaN 
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(c) 



Strain in Ga x ln 1 _ x N grown on GaN 
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FIG. 7. (Color online) Summary of quantities related to or- 
dering in GaaJni-ajN grown on GaN. The SRO parameters 
(a), their difference with respect to bulk case ArQ^ In 
(b), the misfit/relaxation strains calculated on the basis of 
elasticity theory (TOE) or Keating VFF (c), and the criti- 
cal thickness (d). The composition range approximately cor- 
responding to violet-blue-green wavelength of emitted light 
(400-570 nm) is gray shaded. 



the Keating force field is a little bit larger than calculated 
from elasticity. When it comes to the analysis of the 
substrate influence on order parameters, it can be clearly 
seen from Figs. [7] (a) and (b), that the most sensitive to 

(2) 

the presence of substrate is Tq^, as the magnitude of 

(2) 

^-^Gain ^ S * ne largest. Similar behavior is observed not 
only for the GaN, but also for zb-ZnO, and CaO. 

Interestingly, for both negative and positive misfits 

(2) 

(present in the case of zb-ZnO and CaO), the Ar^^ is 
positive, indicating that the presence of the substrate- 
related strain increases the ordering within this shell. 
On the other hand, the quantities Arg^j n and Ar^^ 
are mostly negative, showing the decreasing degree of or- 
dering compared to the bulk case for these coordination 
shells. The magnitude of this effect is, however, lower 

(2) 

than in the case of Ar^^. Finally, the observed behav- 

(3) 

for of ArQ^ In is dependent on the sign of misfit strain. 
For the negative e m i s gt > the decreased ordering is observed 
(Ar^^ < 0), and for positive values of the e m i snt , the 
ordering increases (Ar^^ > 0). 

Let us also mention that the previous calculatio ns by 
Liu and Zunger using epitaxial cluster expansio n 1 31 ^ 32 ! 
predict that already the influence of the GaN substrate 
triggers the phase separation and long-range ordering in 
T = 873 K. They obtained the full miscibility in temper- 
atures above Tm = 1080 K. Our model predicts similar 
phenomena, however, in order to observe them substrates 
related to higher strains are necessary, as discussed in the 
next section. 



B. Large misfits regime 

In this section, we describe the behavior of Ga^In^^N 
on substrates that are related to larger strains, i.e., 
3C-SiC, A1N, and InN. It turns out that in these cases 
the long-range ordering can be triggered. This is detected 
by analyzing the sim-LRO parameter Sq& and its spread 
A^Gai as described in Sec. |IV| For conditions where ei- 
ther jScia or ASca are considerably larger than in the 
bulk case, various ordered structures can be observed in 
the configurations generated during the Monte Carlo run. 
The details of formed patterns depend on composition as 
well as on the magnitude and sign of strain induced by 
the substrate. For the InN base layer, the Ga^Ini-^N al- 
loy undergoes tensile strain, whereas for the 3C-SiC and 
A1N the strain is compressive. Generally, both 3C-SiC 
and A1N correspond to similar misfits as displayed in Ta- 
ble |IV[ therefore, only the 3C-SiC case will be presented 
in detail. The effect of the A1N substrate is analogous. 

Let us begin with the analysis of the InN substrate. 
The complete set of characteristic quantities for this case 
is presented in Fig. [8] Graphs include the dependencies 
on composition of the long-range order parameter £?Ga 
and its spread AS 1 ^, the short-range order parameters 



Gain' 



the relaxation together with misfit strains, and 
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S Ga and A S Ga for Ga x ln 1 _,,N grown on InN 
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Strain in Ga x ln 1 N grown on InN 



0.12 
0.09 



w 0.00 
-0.03 
-0.06 
-0.09 
-0.12 



^relax (TOE) - 


X 












I 


+ 


















+ + 


+ + 












+ + + 


+ + 












































































i 
















I 





(d) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Critical thickness of Ga x ln 1 _ x N grown on InN 



1000 



100 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 8. (Color online) Summary of the quantities related to 
the ordering in Ga^Ini-^N grown on InN. The sim-LRO pa- 
rameter SGa and its spread ASGa (a), the SRO parameters 
(b), the misfit/relaxation strains calculated on the basis 
of the elasticity theory (TOE) or the Keating VFF (c), and 
the critical thickness (d). Open symbols denote the occur- 
rences of the LRO. CH stands for the chalcopyrite and PE for 
the perpendicular planes ordering . The composition range 
approximately corresponding to violet-blue-green wavelength 
of emitted light (400-570 nm) is gray shaded. 




FIG. 9. (Color online) Atomic arrangement in the chalcopy- 
rite structure (the CH pattern). The ideal structure (a) is 
compared with the sample atomic arrangement obtained for 
Gao.5Ino.5N on InN during Monte Carlo simulation (b). Note 
that expansion in the basal (substrate) plane and correspond- 
ing compression in the perpendicular plane is visible, together 
with thermal displacements. 



critical thickness. One can easily see that the sim-LRO 
parameter indicates the presence of the long-range order- 
ing in the range 0.40 < x < 0.90, since either the Sea 
or ASoa reaches high values, compared to the bulk case. 
The observed dependencies of the LRO parameters on 
Monte Carlo time are in these cases similar to the sam- 
ple series presented in Fig.|j(b) and (c). The presence 
of the LRO in the range 0.40 < x < 0.90 is also indicated 
by much higher values of the compared to the bulk 

case. It is worth noting that the whole violet-blue-green 
composition region lies completely in the discussed LRO 
regime for the considered InN substrate. The visual in- 
spection of the structures obtained during simulations 
reveals that for concentrations around x = 0.5 the so 
called chalcopyrite (CH) pattern develops in the exam- 
ined crystal. The ideal chalcopyrite structure is presented 
in Fig. [9] (a) and compared with one of the structures ob- 
tained during simulations in Fig. [9] (b) . Note that even 
though both the deformation due to epitaxial strain and 
the effect of thermal vibrations are visible in the Monte 
Carlo structure, its similarity to the perfect chalcopyrite 
is clear. The concentrations corresponding to the chal- 
copyrite ordering pattern are marked in Fig. [8] with the 
CH label. It is well known from the literature that this 
type of ordering is energetically very favorable in the case 
of strained semiconductor alloys, not only nitrides G22H5H 
Therefore, its appearance in our results is consistent with 
the previous findings. Moreover, the chalcopyrite pattern 
developed in our simulations is oriented perpendicularly 
to the substrate plane (CH_l) as predicted in the recent 
work of Liu and coworkers for large misfits.^ The study 
of the remaining part of the concentration range reveals 
that the CH pattern is not the only ordering option. It 
turns out that different atomic arrangement occurs for 
the higher Ga concentration x > 0.7 . It consists of the 
In planes perpendicular to the substrate. The concentra- 
tions corresponding to this behavior are marked by label 
PE in Fig. [HJ Sample structure with this type of order- 



ing is presented in Fig. 11 (a). Such a behavior can be 



viewed as a certain mode of phase separation. 
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FIG. 10. (Color online) Summary of the quantities related to 
the ordering in Ga^Ini-^N grown on 3C-SiC. The sim-LRO 
parameter Sgh and its spread A£?Ga (a), the SRO param- 
eters FQ aIn (b), the misfit/relaxation strains calculated on 
the basis of the elasticity theory (TOE) or the Keating VFF 
(c), and the critical thickness (d). Open symbols denote the 
occurrences of the LRO of the parallel planes type (PA). The 
composition range approximately corresponding to violet- 
blue-green wavelength of emitted light (400-570 nm) is gray 
shaded. 




FIG. 11. (Color online) Different planar ordering observed in 
Gao.83Ino.17N for two different substrates. For InN substrate, 
the ordering in the direction perpendicular to the substrate 
(the PE pattern) occurs (a). For SiC substrate, ordering in 
the planes parallel to the substrate (the PA pattern) is present 
(b). Only In atoms (minor component) are displayed for clar- 
ity. 



.(4) 

Gain 



This is also indicated by the fact that r^^, Tq^, and 



shift toward high positive values, indicating the 
preference toward Ga-Ga and In-In neighboring. The 
structures without the CH/PE symbols correspond to 
mixed ordering carrying certain features of both arrange- 
ments. The rest of the concentration range, namely for 
x < 0.4 and for x > 0.9, manifests no long range or- 
dering, similarly to the case of GaN, zb-ZnO, and CaO 
substrates. It is also worth noticing that the CH ordering 
is observed close to the limit of elastic accommodation of 
strain regime, marked by the existence of critical thick- 
ness solutions in Eq. (111. The PE ordering in turn is 



observed well beyond this regime. Therefore, in reality 
the PE case might be difficult to observe due to the very 
poor crystalline layer quality in this misfit region. 

The behavior of Ga x Ini_ x N on substrates correspond- 
ing to large compressive strains also leads to the forma- 
tion of the long-range ordering. The detailed dependen- 
cies on concentration for the long-range order character- 
istics Sea an d ASca, the short-range order parameters 



-.(0 

Gain' 



the relaxation and misfit strains accompanied by 
the critical thickness for SiC substrate are presented in 
Fig. [TUj The analogous results for A1N are not included, 
since they were almost identical. From our simulations 
it emerges that for the high gallium concentrations in 
the range 0.60 < x < 0.85 the long-range ordering can 
be observed. Even though Sg& remains close to zero, 
the dependence of ASca on concentration experiences 
abrupt change in the ordering region. This change indi- 
cates that the Sea dependence on MC time is similar to 
the sample presented in Fig. 3] (c) . This is also confirmed 
by the larger values of the short-range order parameters 
in this regime. It is worth noting that all become 
positive, indicating Ga-Ga and In-In neighboring prefer- 
ence. Visual inspection of configuration from the region 
with LRO reveals that indeed atoms tend to cluster in 
the planes parallel to the substrate. This is depicted in 



Fig. 10 as the PA pattern. Such a behavior can be inter- 
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preted as a certain mode of phase separation. The sample 
structure corresponding to this arrangement is displayed 
in Fig. 11 (b). For the 3C-SiC substrate, this ordering 
occurs partially within the region of elastic strain accom- 
modation marked by the existence of solutions for the 
critical thickness model from Eq. (11), so it should be 



possible to access it experimentally. Moreover, the whole 
discussed area of the LRO existence lies within the tech- 
nologically important violet-blue-green composition re- 
gion. For the Ga concentrations below the x = 0.6, which 
correspond to the very high compressive strains, no sig- 
nificant indications of the long-range order are observed. 
Interestingly the parameters T^ In , T^ In , and T^ In in 
this concentration region have absolute values lower than 

(2) 

in the unstrained bulk case. The Ig^, in turn, simi- 
larly to the moderate strain case, deviates mostly from 
the bulk ordering and its absolute value increases. These 
phenomena, however, occur well outside the elastic ac- 
commodation of strain regime, where in reality extremely 
poor crystal quality would be obtained. 



VII. SUMMARY 

In this paper, the ordering phenomena in ternary ni- 
tride alloys have been examined in a great detail. We 
have investigated how temperature, biaxial strain related 
to the presence of substrate, and the change of compo- 
sition in the whole available range could influence type 
and degree of ordering. Our approach included vibra- 
tional contribution, which is not so common in the lit- 
erature, where restricting to configurational degrees of 
freedom only is a frequent approximation. As a start- 
ing point bulk Ga^Ini_ x N, ALTni-^N, and Al^Gai-^N 
have been investigated. The ALGai-^N has been found 
to follow closely uncorrelated random alloy picture, with- 
out both SRO and LRO. For the mixtures containing 
indium, no signs of long-range ordering or precipitation 
have been found in the bulk case; however, a consider- 
able degree of short-range ordering has been observed. 
The short-range order parameters corresponding to the 
first four coordination shells followed the — / + / + /— 
sign pattern, which agrees with the behavior observed in 
simulations for Ga^Ini-^N on the basis of cluster expan- 
sion modelPS Both materials deviate the strongest from 



the random uncorrelated alloy at high indium concen- 
trations around 65%-70%. This means that theoretical 
methods neglecting SRO should yield the largest system- 
atic error in this concentration range. When it comes 
to resolvi ng the I n clustering controversy reported in the 
literature) 5 * 7 * 9 * 1 ^ our findings could provide a certain sup- 
port for the experimental results showing uniform dis- 
tribution of In atoms in Ga^Ini-^N samples. However, 
the model of bulk alloy might not be fully appropriate 
here, since examined nitride alloys are usually employed 
in (opto)electronic devices in the form of strained epi- 
taxial layers. Therefore, the next phase of our studies 
has been devoted to the influence of the substrate re- 
lated biaxial strain on Ga^Ini-^N layers. A number of 
different substrates has been examined and classified into 
two groups. The first group, associated with moderate 
strains, contained GaN, CaO, and zb-ZnO. In this case, 
in the examined temperature T = 873 K, still only SRO 
has occurred. The second and more interesting group 
comprised substrates yielding larger misfits, namely A1N, 
3C-SiC, and InN. For these substrates, the appearance of 
the LRO has been detected by introduced sim-LRO pa- 
rameter. For compressive misfit strain (A1N, 3C-SiC), 
the parallel planes (PA) pattern has developed for the 
In content 15%-40%. For tensile misfit strain (InN) the 
LRO regime has covered the region of 10%-60% of In con- 
tent. The chalcopyrite ordering (CH) has emerged, when 
the indium concentration was around 50%, whereas for 
lower values down to 10%, the perpendicular planes (PE) 
pattern has been found. The PA and the PE cases rep- 
resent certain types of precipitation. The above results 
shed more light on the issue of ordering in nitrides oc- 
curring in thermodynamical equilibrium. We believe that 
the knowledge presented here can facilitate further mod- 
eling, e.g., of electronic structure in nitride alloys. It can 
be also be helpful in the interpretation of experimental 
findings. 
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